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SECTION  I 


THE  MODEL 

In  estimating  the  performance  of  an  error-correcting  device  on 
a  specific  communication  channel,  it  is  necessary  to  find  a  meaning¬ 
ful,  yet  tractable,  mathematical  model  for  that  channel.  An  examination 
of  data  from  real  channels  suggests  that  most  channels  pass  through 
three  distinct  phases.  The  first  phase,  which  nearly  any  error- 
correcting  scheme  can  handle  successfully,  is  that  of  long  periods  of 
practically  error-free  transmission.  The  second  phase,  which  is  the 
antithesis  of  the  first  phase  in  that  no  scheme  can  expect  to  correct 
it,  is  that  of  complete  loss  of  signal  for  substantial  periods  of  time. 
The  third  phase  might  be  described  as  a  generalized  sputter  or  bursts 
of  error  bursts,  and  it  is  this  phase,  if  it  occurs  frequently,  for 
which  error  correctors  should  be  designed.  This  third  phase  may 
sometimes  be  described  by  means  of  a  two- state  model.  The  three 
phases  suggest  a  Markov  process  with  four  states,  but  such  a  process 
is  mathematically  unwiedly.  However,  when  the  three  phase  picture  is 
reasonably  correct  and  the  sputter  phase  is  accurately  modeled  for 
purposes  of  estimating  coder  performance,  it  is  a  simple  matter  to 
make  corrections  in  such  estimates  for  those  periods  when  transmission 
is  nearly  error  free  or  when  the  signal  is  lost. 

The  use  of  such  a  two- state  model  was  suggested  to  us  by  some  work 
by  Gilbert  [1].  Gilbert  describes  a  model  for  binary  error  distributions 
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in  channels  subject  to  noise  bursts.  Let  {x.}  be  the  error  process  with 

x#  =  1  for  an  error  in  the  i^  demodulated  bit,  x.  =  0  for  no  error, 
i  l 

Two  states,  designated  G  and  B,  of  the  channel  are  postulated  such  that 

th 

at  the  i  bit  the  state  S.  is  G  or  B,  and  the  state  S4|,  at  bit  i  +  1 

l  l+l 

depends  only  on  Thus,  the  state  sequence  {S^ (  is  a  simple  Markov 

chain  described  by  the  two  transition  probabilities,  P:  G  ->  B  and 
p  B  -»  G.  We  also  use  Gilbert's  notation  Q  =  1  -  P  and  q  =  1  -  p  for 
the  G  -*  G  and  B  -*  B  transitions  respectively  Let  h,k  denote  respectively 
the  probabilities  of  a  correct  bit  in  B  and  in  G  Then  P[x^=l (  =  1  -  h  -  h 1 
if  S.  -  B  and  P{x.=l]  -  1  -  k  -  k*  if  S.  -  G 

l  l  i 

Using  this  model  with  k  =  1,  Gilbert  obtained  a  good  fit  to  certain 

phone  line  error  data.  The  statistic  fitted  was  the  probability  of 

occurrence  of  zero  (=  error  free)  runs  of  length  at  least  K  (The  use 

of  k  <  1  has  also  been  considered  by  Elliott  [5].) 

At  MITRE,  we  have  found  sets  of  values  for  the  model  parameters 

P,  p,  h,  k,  which  yield  good  fits  to  error  data  for  several  different 

types  of  communication  media  One  important  statistic  involved  is  the 

probability  of  specific  error  densities  in  various  length  blocks. 

Given  that  we  have  p,  P,  h,k,  [2,3,4],  we  present  the  recursive 

method  used  to  calculate  the  probability  of  K  symbol  errors  in  a  block 

of  n  symbols  where  each  symbol  consists  of  m  bits. 
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SECTION  II 


RECURSIONS  -  PART  I 

In  Appendix  II  of  [2]  (or  Appendix  B  of  [3]  or  [4])  we  presented  a 
brief  outline  of  the  recursive  method  used  to  calculate  the  probability 
of  K  symbol  errors  in  a  block  of  n  symbols  where  each  symbol  consists 
of  m  bits  Now  we  present  the  outline  in  full.  Since  the  last 
documentation,  the  technique  has  been  extended  to  permit  the  n-symbol 
blocks  to  be  interleaved  or  time  spread.  (To  interleave  s  blocks  means 
to  transmit  sequentially  the  first  symbol  of  each  of  s  blocks  followed  by 
the  second  symbol  of  each  of  those  blocks,  etc.  Thus,  a  given  m-bit 
symbol  is  transmitted  s  symbol  times  after  the  symbol  preceding  it  in 
its  block . ) 

Letting  T  and  U  represent  either  of  the  states  G  and  B,  we  define 
TOU(t)  =  P(xx  ■  •  •  •  ■  ■  0  and  S^_  =  U|SQ  =  T) 

TlU(t)  =  P(for  some  i  £  t,  x  ■  1  and  c  U|Sq  =  T) 

Then 


G0G(1)  =  Qk 

G0B(1)  =  Ph 

G1G(1)  =  Qk' 

G1B(1)  =  Ph' 

B0B(1)  =  qh 

B0G(1)  =  pk 

B1B(1)  =  qh' 

B1G(1)  =  pk' 

and 

GOG(t)  =  [GOB (t- 1 )  p  +  GOG(t-l)  Q]  k 

GlG(t)  =  [GOB(t-l)  p  +  GOG(t-l)  Q]  k'  +  GlB(t-l)  p  +  GlG(t-l)  Q 
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GOB(C)  =  [GOB(t-l)  q  +  GOG(t-l)  P]  h 

GlB(t)  =  [GOB (t- 1 )  q  +  GOG(t-l)  P]  h'  +  GlB(t-i;  q  +  GlG(t-l)  P 

BOB(t)  =  [BOB(t-l)  q  +  BOG(t-l)  P]  h 

BlB(t)  =  [BOB(t-l)  q  +  BOG(t-l)  P]  h'  +  BlB(t-l)  q  +  BlG(t-l)  P 

BOG(t)  =  [BOG(t-l)  Q  +  BOB(t-l)  p]  k 

B1G(C)  =  [BOG(t-l)  Q  +  BOB(t-l)  p]  k'  +  BlB(t-l)  p  +  BlG(t-l)  Q 

Shortly  after  the  program  was  written,  we  discovered  that  GOG(m), 
GOB(m),  GlG(m) ,  GlB(m),  BOG (m) ,  BOB(m) ,  BlG(m) ,  BlB(m)  could  be  obtained 
from  a  difference  equation  in  powers  of  J  and  L  (see  below).  Since  on 
our  computer  (IBM  7030),  it  took  under  a  second  to  compute  all  eight 
quantities,  we  decided  not  to  use  the  difference  equation.  However,  we 
work  out  two,  and  give  all  eight  results. 

GOG(t)  =  [GOG(t-l)  Q  +  GOB (t- 1 )  p}  k 
GOB (t )  =  [GOB(t-l)  q  +G0G(t-l)  P}  h 

The  eigenvalues  come  from  the  2nd  order  linear  difference  equation: 

ft+l  ‘  (Qk  +  qh)  ft  '  (P'Q)  ft-l  =  °* 
that  is,  2j  =  Qk  +  qh  W(Qk  +  qh)2  +  4hk(p-Q) 

and  21.  =  Qk  +  qh  -  ^(Qk  +  qh)^  +  4hk(p-Q) 

Thus  we  have 

GOG(t)  -  ^  +  0^  Lfc  , 

GOB(t)  =  fil  l'1 
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and  we  get 


<*!>  0 '2,  ^  and  ^ 


from  the  initial  conditions 

GOG(O)  =  1,  G0G(1)  =  Qk 
and  GOB(O)  =  0,  G0B(1)  =  Ph  . 

So  al  +  °2  =  +  O^L  =  Qk,  which  yields 

GOG(t)  =  {(Qk  -  L)/(J  -  L)}  Jl  +  {(J  -  Qk)/ (J  -  L)}  L* 

Also,  +  ^2  =  ^1J  +  ^2L  =  Ph>  which  yields 

GOB(t)  =  {Ph/(J  -  L)}  (JC  -  LC) 

All  eight  solutions  are  as  follows: 

GOG (m)  =  {(Qk  -  L)/(J  -  L)}  J™  +  { (J  -  Qk)/(J  -  L)}  L™ 

GOB(tn)  =  {Ph/(J  -  L)}  (J™  -  L™) 

BOG (m)  =  [pk/(J  -  L)}  (J™  -  L™) 

BOB(tn)  =  {(qh  -  L)/(J  -  L)}  J™  +  { (J  -  qh)/(J  -  L)}  L™ 

GlG(tn)  =  p/(p  +  P)  +  [P/(p  +  P)}  (Q~p)m  -  {(J  -  qh)/(J  -  L)}  J1" 
-  {(qh  -  L)/(J  -  L)}  Lm 

GlB(tn)  -  ^  [1  -  (Q~ P )m]  -  fl I  ^ 

BlB(m)  =  —j  +  (Q-p)m  -  ^  +  !Qk-LjLm] 

BlG(tn)  =  [1  -  (Q-P)™]  -  ^  (J1"  -  L™) 
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SECTION  III 


RECURSIONS  -  PART  2 


Again  letting  T  and  U  represent  either  of  the  states  G  and  B, 


we  define 


x 


-  0,S 


sm 


T) 


=  P(m-bit  symbol  after  s  symbol  times  is  correct  and 

ends  in  state  ujstate  T) 


TlUI(s)  =  P(for  some  1  £  i  £  m,  x/  ...=1,S  =  u|Sn  =  T) 

(s-I)mfi  sm  1  0 


=  P(m-bit  symbol  after  s  symbol  times  has  at  least  one  bit 

error  and  ends  in  state  u|state  T) 


Let  GXG  =  GOG(m)  +  GlG(m) ,  GXB  -  GOB (m)  +  GlB(m),  BXG  =  BOG(m)  +  BlG(m), 
and  BXB  =  BOB (m)  +  BlB(m) 

Then  GOGI(s)  =  GXG  GOGI(s-l)  +  GXB  *  BOGI(s-l)  and  BOGI(s)  = 

BXG  GOGI(s-l)  +  BXB  *  BOGI(s-l). 

(There  will  be  similar  equations  in  G1BI,  B1BI  and  GOBI,  BOBI ,  and 

G1BI,  B1B1,  but  they  have  the  same  eigenvalues  and  will  be  omitted.) 

TOGI(s)  -  [GXG  +  BXB]  TOGI(s-l)  +  [GXG* BXB -GXB  BXG]  T0GI(s-2)  «  0 

This  yields  the  eigenvalues: 

2  CT  =  GXG  +  BXB  +  >/[GXG  +  BXB]2-4  [GXG  •  BXB  -  GXB  ■  BXG] 

2  t  =  GXG  +  BXB  - +  BXB]2-4  [GXG  '  BXB  -  GXB  •  BXG] 
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Since  GXG  +  GXB  =  BXB  +  BXG  =  1, 


a  T  =  GXG  •  BXB  -  GXB  *  BXG  =  (1-GXB)  (1-BXG)  -  GXB  •  BXG  =  1-GXB-BXG 

a  +  T  =  GXG  +  BXB  =  2  -  GXB  -  BXG  =  1  +  a  T  . 

Hence  a  =  1 ,  and  T  =  1  -  GXB  -  BXG  =  GXG  +  BXB  -  1 . 

Thus  we  have 

GOGI(s)  =  X1  1S  +  \2tS 
and  BOGI(s)  =  1S  + 

and  we  get  and  from  the  initial  conditions. 

G0GI(1)  =  GOG(m),  G0GI(2)  =  GXG  *  GOG(m)  +  GXB  •  BOG(m) 

BOGI(l)  =  BOG(m),  B0GI(2)  =  BXG  *  GOG(m)  +  BXB  •  BOG(m) 

Hence 

A  1  +  X2  '  T  =  G0G(m^>  \  ’  l2  +  ’  T2  =  GXG  *  G0G(m)  +  GXB  *  BOG(m). 

Solving,  X  =  [ GOG (m)  (1  -  BXB)  +  BOG (m)  •  GXB]  /  (1  -  t) 

\2  =  [GOG (m)  BXB  -  BOG(m)  GXG]  /  (1  -  t)  =  [GOG(m)  -  ]  /  t 

Also , 

H  •  1  +  ^2  ‘  T  =  BOG(tn),  ^  •  l2  +  ^  •  T2  =  BXG  *  GOG (m)  +  BXB  •  BOG(m) . 

Solving,  =  [BOG (m)  (1  -  GXG)  +  BXG  *  GOG(m)]  /  (1  -  t) 

42  =  [BOG(m)  -  /  T 


7 


Therefore,  GOGI(s)  =  X^  '  Is  +  X^  *  tS  ,  and  BOGI(s)  =  Mj'lS  +  ti2  *  TS , 
where  X^ ,  are  given  above.  Since  s  is  fixed  for  any  given 

application,  we  will  refer  to  GOGI(s)  as  GOGI ,  and  to  BOGl(s)  as  BOGI. 

Consider  the  first  i  m-bit  symbols  of  a  random  interleaved  block 
Let  GB(i,j)  =  P(j  symbol  errors  in  i  symbols  and  b|Sq  *  G). 

Similarly,  we  define  GG(i,j),  BB(i,j)  and  BG(i,j). 

Then 


GG(1 ,0)  -  GOGI 

GB(1,0)  -  GOBI 

GG(1 , 1 )  =  G1GI 

GB(1 , 1 )  =  G1BI 

BG ( 1 , 0 )  =  BOGI 

BB (1 ,0 )  =  BOBI 

BG(1 , 1 )  =  BlGI 

BB  (1 , 1 )  =  BlBI 

Finally,  for  i  =  2,  ■  ,  n  and  j  =0,  '  '  ,  i 

GG  (  i  ,  j  )  =  GG  (  i  - 1  ,  j  .)  GG  (1,0)  +  GB(i-l,j)  BG(1,0; 

+  GG(i-l.j-l)  GG(1 ,1 )  +  GB (i- 1 , j -1 )  BG(1,1) 

GB ( i  ,  j )  =  GG ( i - 1 , j )  GB (1 , 0 )  +  GB(i-l,j)  BB(1,0J 

+  GG ( i - 1 , j - 1 )  GB (1 , 1 )  +  GB (i- 1 , j - 1 )  BB(1,1) 


BG(i  ,  j )  =  BG(i-l , j )  GG (1,0)  +  BB(i-l,j)  BG(1,0) 

+  BG(i-l , j-1)  GG (1,1)  +  BB (i- 1 , j - 1 )  BG(1,1) 

BB (i  ,  j )  =  BG(i- 1 , j )  GB ( 1 , 0 )  +  BB(i-l,j)  BB(1,0) 

+  BG(i-l , j-1)  GB (1  , 1 )  +  BB(i-l  ,  j )  BB(1,1). 
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Finally , 


P (random  bit  is  in  G)  =  a.  =  and  hence 


P(K  symbol  errors  in  n  symbols  with  s  blocks  interleaved) 


CL  [GG(n,K)  +  GB(n,K)]  +  (1  -  a)  [BG(n ,K)  +  BB(n,K)] 
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SECTION  IV 


INPUT  FOR  THE  PROGRAM 


XM 

XN 

XNEST 


XIPER 

IK 

CP 

SP 

H 

SK 


No  of  bits/symbol 
No.  of  symbols/block 

Largest  number  of  symbol  errors  to  be  considered 
(if  this  field  is  blank,  XNEST  =  XN) 

Number  of  interleaved  symbols  (0  or  1  means  1) 

0  or  not  equal  to  0  (0  means  continue  with  CP,  SP,  H,  SK; 
not  equal  to  0  means  read  new  parameters) 

P 

P 

h 

k 


» 


« 
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SECTION  V 


OUTPUT  OF  THE  PROGRAM 

TIMEX  =  A8  representation  of  the  time  read  from  IBM  7030  Time  Clock 
by  a  STRAP  coded  routine  (one  can  call  his  routine  or  omit 
it  altogether) 

CP,  SP,  H,  SK  as  above 

ALPHA  =  P  (random  bit  is  in  G) 

M  (or  XM)  as  above 

N  (or  XN)  as  above 

NS  =  No.  of  recursion  terms  to  be  attempted 
WPI  =  P  (symbol  error) 

WY155  =  P  (no  errors  when  s  =  1) 

WPMU2  =  mean  number  of  errors  in  a  block 
FCMEAN  =  mean  number  of  errors  given  an  error  occurred  (when  s  =  1) 
For  J,  L,  A,  B  see  [2,3,^] 

IS  =  number  of  errors 
P  *  P  (IS  symbol  errors  in  a  block) 

R  =  P  (IS  symbol  errors  in  a  block |error  occurred) 

Q  *  P  (£  IS  symbol  errors  in  a  block) 

QH  *  P  (>  IS  symbol  errors  in  a  block) 

S  *  mean  number  of  bits  between  blocks  with  >  IS  symbol  errors 

PBAR  *  contribution  to  mean  number  of  errors  per  block  made  by 

probabilities  actually  calculated  (if  one  wants  the  whole 
mean,  then  NEST  =  0  or  N;  the  same  applies  to  VAR,  SVAR, 
and  CMEAN) 
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VAR 

SVAR 

CMEAN 


contribution  to  variance  about  PBAR  made  by  probabilities 
actually  calculated  by  the  recursion 

=  approximate  standard  deviation 

=  contribution  to  mean  number  of  errors  given  an  error 
occurred  in  the  block  made  by  probabilities  actually 
calculated  by  recursion 
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APPENDIX 


PROGRAM 


V  fcftSIdH  5572R/G5  -  CORRECTION  LE¥Etl  £T3/2B?i5 


00000 

0.QU0I 


C...T0  CALCULATE  the  recursion  probabili ties 
COMMUN  771/  TIMEX 

DIMENSION  GGI515) t  G6{ 515) ,  BGj515_L»  OBC515) 


00005 

00006 


_ _ _ _  ■  -^Sfll2Crtl  ,  G1BI2031  = 

XHl50I,Wl5ai,XNESI  (501*  n=T^XTF*iU  50) 

404  REA01701,  LK 

1701  FORMAT  ill) 

C...KEAO  IN  THE  NO.  OF  BITS  PER  SYMBOL,  THE  NO.  OF  SYMBOLS, 


jswmr~  7gr 


00009 

00010. 


C. .  .R 

604 

702.. 


HO*  MANY  SYMBOLS  TO  GO  THRU  *  HOk  MANY 
Ll-tr 

TEQgMAT  IIMX,  16,  2Xt  13,  .  _  .  .141= 

EAD  IN  CP,  SP,  H,  SK 
READ  702,  IK,  CP,  SP,  H,  SK 
FORMAT  (11,  E17.10,  3CE18.il)) 


iTitKi 


m^z 


m  Z&&  K*l,  LK 


M  XMTK  T 

|  XNIK) 


00014 

00016 

00016 


NEST  -  XNEST(K) 
IPERD  =X  IPER ( K ) 
IF(NEST.EQ.O)  NEST 


000  IT 
0B31R 


,DQSTS 


NS  *  NEST  +  1 

CALL  TIME 

PRINT  300,  TIMEX 


00020 


00021 


300 


00025 

00026 

00027 


FORMAT ( 1H 1 , Ad  ,21X,  77H*CALCULAT1UN  OF  GILBERT  CHANNEL  ERROR  PROBAB 
llLlTItS  USING  RECURSION  FORMULAS*/////) 

MUD  _= _ M*N_ _ _ _ _ _ 

fMUO  -  FLOAT  (MUD) 

MU01  HUD  -  1 

_ =^  !>  -  H  _  1 _  .  _  .  _ _ 

SKCUMP  =  1*  -  SK 

CQ  =  1  *  -  CP 


so 


1  •  -  SP 


GAMMA  =  1./1CP4-SP) 
ALPHA  »  SP*GAMMA 
BETA  =  CP*GAMMA 


00031 

00032 

00033 


301 


PI  =  ( SP*  SKCUMP  ♦  CP*HCOMP) *GAMMA 
PRINT  301,  CP,  SP,  H,  SK,  ALPHA 
FORMAT  ( 5X»  2  7H0HANNE L  PARAMETERS 


CP  =  , E 1 1 • 6 ,  5X ,  5HSP  * 


4ElU64Xf  H  *  f  F8.6,  5X ,  5H  K  =  ,Elfc*lG*S3U  mALPHA  =  *  FBJA) 

PRINT  302*  Mt  N*  PiSi  Pi*  lPEfUl 

ZfitfflE  F AR&fgTSRS-  =  M-=  +-ZW+  -131IT  W  " 


1 1 2 X ,  35HN0.  OF  RECURSION  TERMS  REQUESTED  =  , 

2  Eli. 7,  75X ,  14, 

3  $,/) 


15/  7X ,  SP1  =  S, 

$  BLOCKS  INTERLEAVED 
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* 


0CC36 

00037 

00038 

mN 

WT 

HR 

* 

SP  -  CQ 

8K*CQ  +  H*  SO 

WT* *2  €  4.*H*SK*WN 

f  , 

00039 

00040 

MSUKT 
w  J  = 

=  SORT (HR) 

t  * T  ♦  HSOR  T ) /2 . 

00041 

00U42 

WL 

WA 

(WT  -  wS0RT)/2. 

WJ  ♦  wN  *(CP*Sk*HCURP  +  SP*H*SKCOMP  )  / 

1  (CP*HCUMP  +  *P*$KCUMP  ) 


UQC43  WA  =  HA  /  WSQRT 


00044 

W  8  = 

HA  -  1  . 

00u45 

H  X  2  = 

( WA* ( 1*  — WJ**M  J /( 1 .  — W  J  J )  -  ( w  B*  (  l.-WL**M)/( l.-WL )  ) 

0  0  0  4  6 

00G4  7 
00048 

HP  I  = 
HPMU2 

K  J  * 

PI  *  HX2 
=  FLOAT (N )  *  HP 1 

F  LG AT  t  MUO )  *  ALGGlO(WJ) 

0UG49  KL  -  FLUAT(MUD)  *  ALUGlOIwU 


00050 

1FIR0. 

LT. 

-50.  ) 

F 

=  0* 

00051 

IFIRJ. 

GE. 

-50. ) 

F 

=  WJ**MUO 

000  52 

(  F  ( R  L  . 

LT. 

-60.  ) 

G 

=  0. 

00053 

(Ft*L  . 

GE. 

-50.  ) 

G 

=  WL**MUL> 

00054 _  HY155  =  1.  -(  trtA*  {  1  .-F  )  /  t  1  .-HJii-HdM  1  .-G )  /  1  U-HL  )  )  *P  1 

00055  fCrtfcAW  =  WPMU2  /  (1.  -  WY155) 

00056  PRffSiT  806,  WPJ,  W  Y  1  5  5  »  UPHU2  , FbMfc  AN ,  W  J  ,  *L*  WA,  Wtt 

00057 _  806  FORMAT  (IX,  20HPRQ6  SYMBOL  ERROR  =  ,£12,6*  1H*,  2X*  I7HPRQB  jNO  £R 

1K0RS  =  ,  b  1  2 . 6  *  1 H*,  2  X ,  OH  Mt  AN  =  ,tl2.6,  1H*,  2X,  SCUNOITlUNAL  ME 

2AN  =  i,  1 1  2  •  6  ,  $*$  /  1JX,  4HJ  =  ,  3  X ,  £16.10,  1H*, 

3 1  5  X ,  <♦  HL  =  ,  3X ,  F9.6,  1H*  ,  6X ,  4H  A  >  ,  3X,  F9.6,  lH*,  20X,  4HB  « 

4fF9.6t  1H*///1 


00058 

GC059 

G0C6C 

00061 

00062 

00063 

00064 

00065 

00066 

00067 

00068 

0006* 
00070 
00071 
00072 
00073 
00074 
£00075 
00076 
00077 
00078 
0  00  79 
00080 
00081 


GOG( 1) 
GlG(l) 
GOB ( 1) 
Glbl 11 
BJBl  1 J 
BIB ( 1 > 
BOG ( 1) 
B1G(1) 

I  Ft  M.EU 
00  50 
J 1  =  J- 


=  CO*  SK 
*  €0  - 
=  C  P*  H 
=  CP  - 
=  SQ*H 
=  SQ  - 
=  SP*  SK, 

=  SP  - 
li  GO  TU 
» M 


G0G(1> 


GOB ( 1 J 
BOB ( 1 ) 

BOG (1) 

8u 


50 

80 


TA 

GOGtOJ  * 
GIG ( J )  = 
TB 

GOBtJ)  = 
G1B(J)  « 
T£ 

BOB ( J )  = 
B1BU1  = 
TO 

BOGtJ)  = 
B1G(J)  = 
PRINT  3C3 


*  G1G(JU*CQ 
G  OG  t  J 1 )  *  CP 


*  Glri(Jil*5P 


J  =  2  , 

1 _  ... _ _ _ 

=  G0b("J1)*SP  +  G0GIJ1)*C0 
»  T  A*  SK 
=  TA  -  GOG(J) 

=  G OB t  J 1 )*S0  < 

=  T  B*H 

=  TB  -  GOB ( J ) 

=  BQd(Jl)*SO  ♦  BOG ( Jl) *CP 
=  TC*H 

TC  -  BOB ( J ) 

BOG( Jl >*CW 
TO*  SK 

TO  -  BOG(J)  ♦  BIG  (  J 1  >  *C0  +  BlbtJi  )  *S  P 
GOG(M),  GlG(H),  GOB (Ml  *  GlB(rt) 


♦  GIBt Jl) *SQ  +  G1G( Jl ) *CP 


*  BIB ( 31 ) *S0 
B  OB (  Jl J  *  SP 
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00092  303  fOKMAt  6HG0G  ?  ,  bl2«5*  5X?  6Hult»  =  *  £I2»5,  5X^  6MG08  —  y  - 

=-^=-=^  SXg^ttOig_-i«>  ♦  €1 2, 3*  _  -  - - ~-^== 

000B3  PRINT  304,  BOG ( M )  ,  BIG(M),  B08(M),  BIB ( M) 

000B4  304  FORMAT  (5X,  6HB0G  =  ,  E12.5,  5Xf  6HB1G  =  ,  E12.5,  5X,  6HB0B  *  * 

IE  12 • 5 ,  5X  f  6HBIB  =  ,  E12.5//I 

00085^^^  ^TRP I 1  -  ALPHA*  (GlSTMf  *  G1G1H)  )  *  dETA»(  BIB  i  *i  +  0101*0 

00086  .  PRINT  505,  RPII 

QOOB7—  —  505  fO^HAf  JiX*  30HPROB  RANDOM  SVM80L  IN  ERROR  ■  ,  EL2«6,  1M+/// } 

OOOBB  CALL  TIME 

000B9  GG( 1)  =  GOG(M) 

00090  GG( 2 )  =  G1G(M) 

80091  GBC1I  -  G08f«l=^=7== 

jSg6«  —  g*t2>  - 

00094  BG( 2 )  =  BIG(M) 

00095  BB<  1 )  =  BOB  (  M) 

00096  BB ( 2 )  *  B1B(M) 

00092-  IF  ( fPEKD  ,N£*  01  GO  TO  9001 

00098  GOG  I  =  bOG  f  M I 

00099  GlCI  *  GIGIMI 

00100  GOBI  =  GOB (M ) 

00101  G IB  I  =  GIB(M) 

00102 _ BOGI  =  BOG(M) _ _  _ 

00103  BIG!  =  BiGIM) 


00107  9001  GXG  =  bGU)  ♦  GGt  2 1 
00108  GXB  ^  GB ( I >  4  Gb ( 2  I 
00109 _ 6XG  =  6GI 1 )  *  BGI  21 


00119  PRINT  9003,  SIG,  TAU 
00120  9003  FORMA!  ( 6X ,  SSIGMA  = 
-00121  AG *  (  GG<  1 )  *GT * 


$,  E20.12,  10X,  STAU 
GXB  *  BG(1))  /  DEN 


$,  E20.12/' 


mx  ar 

-«*» 

80124 

00125 

00126 

00127 


GOGI  =  AG 
BOGI  =  AB 
AG  =  ( GG( 2 ) 


*GT 


CG 

CB 

♦ 


♦  G 

♦  G 
GXB 


♦  BG  (  2) )  /  DEN 


00129 


AB 


fGGf21  -  AG  I  /  TAU 

10X0  *  00£2t  *  BGI2I  *  8T)  /  0£H 


% 
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0G130 

00131 

00132 

0  C  1  33 

0013* 

00135 

00136_ 

0C137 

0013B 

00139 

001*0 

0C141 

00142 

001*3 

00144 

001*5 


CB  =  IBg(2)  -  AB)  /  TAU 
G1GI  =  Ab  ♦  CG  *  G 

ttlGI  =  AB  ♦  CB  *  G 

AG _ =  <_GB1 1)  *GT  ♦  GXb 

.CG  =  (GBill-AG)  /  TAU 
AB  =  (BXC  »  Will  +  BttUl 
CB  =  IBBi 1)  -  AB)  /  TAU 


*  B6< 1)1  /  OEN 


*  BT)  /  O £H 


OEN 


CG  *  G 
C  B  *  G 

*6T  _♦ _ GXB  ♦  BB ( 2 ) )  / 

CG  =  ( G8(  2)  -  AG)  /  TAU 

AB  *  <8XG  *  Gbl2)  ♦  Bb(2)  *  BU  /  D£N 

CB  =  IBB12)  -  AB)  /  TAU 


GCB I  =  Ab 
BOBI  =  AB 
AG  =  ( gB ( 2 ) 


00146 


00161 


G1BI  =  Ag  ♦  Co  *  G 

B1BI  =  AB  ♦  CB  *  G 

_ PRINT  9003  t  GOG 1 t  G1G1 >  GOBI  t  GlBl t  BOG  1 »  B1GI ,  BOB  It  B1B1 _ 

900!*  £GitMAT  I/sX*  7MG0GI  =  WE12^5*  4X*  7HG1GI  ,  £12.5*  4X,  7HGU&I 

i  €12. 5«  *X»  7HG10I  *  t  £12.5/  5Xf  7HB0G1  •  ,  £12.5*  *Xf  7hB1gI 

7H30B I  =  ,  £12.5,  *X,  7HBIS1  =  ,  E12.5//1 


Z  £12.5*  4X i 


p  -*! 


00147 

C 

00146 

9002  1 F ( NEST . EQ . 

. .  .CALCULATE  THE 
DO  51  L 1 =2 

00149 

LZ  = 

MING 

00150 

LY  = 

Cl 

00151 

IFIL1 

.GT.  1 

00152 

IFILI 

.6T.  1 

00153 

N 1  = 

LI  ♦  1 

00154 

GG(N1) 

»  N£$T  +  1 


Go  I L  1 ) *6 1G I 


GB ( LI ) *B1G ( 


00155  G&IN1) 

00156  BGiNU 

00157  BBlril)  • 

00158  97  CONTINUk 

00159  DU  02  L2=2,LY 

00160  N2  =  LY  -  L 2 


GG (LI )  *61B I 
BG  <  L 1 1 *u 1g I 
BG I L 1 ) *6  IB  I 


GSTL1)*B1BI 

MTiLU^GI 

BBt€n*BlBl 


+  2 


N3  =  N2  -  1 


00162  YAH1  *GGIN2)*G0Gi 

00163  GB4N2 ) =6B I NZ) *BOBI 

00164  G&TN2)  =  YAH1 

00165  YAH2  =BG(N2)*G0GI 

00166  BB(N2)=BB( N2)*B0BI 

00167   52  BG ( N2 )  =  Y  AH  2 


+GGIN3)*G1GI 
♦CB (N3T*diBl 


+GSTiiZT4*  O&i 


+&6iN3}*B10I 

♦GGW>*GlBi 


00166 

00169 

0C17C 


♦Bo 1N3 ) *6161 
♦BB(N3)*BlbI 


♦  BB  ( N2 ) *bOg I 
♦ BG 1 N2 ) *G0B I 


♦  BB ( N3 ) *B1G I 

♦  BG ( N  3 ) *G1B I 


YAH3 


zSftlll  * 

&&1I1  — *= 


=££tU*GQGX 
6B( 1 1  *B08 i 
YAH  3 


♦  bb(  1 )  *  BOG  I 
♦  Gb(ll  *0081 


0C171 

00172 

00173 


YAH* 
BBC  1) 
BG(  1 ) 


BG  (  1  I 
=  BUll) 
*  YAH* 


♦GOG  I 
♦BOB  J 


BBU) 
BG  ( 1 ) 


*B0GI 
♦  GOBI 


51  CONTINUE 
90  CONTINUE 


{FINEST. 01 .33)  PRINT  503 
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00177  5  03  FORMAT  (  /  //  //  //  //////  / 50  X*  33HI  SYHBUL  DA  T  A  STARTS  UN  NEXT  PAGE  )  ) 

00178  IFtNbST  •  GT .35)  PRINT  502 

00179  SO 2  FORMAT  ( 1H1) 

06180  '  T»klNl  506,  TIMEX 

00181  500  FORMAT I  lx , A8/  IX,  SSYMBUL  NUMBtftSlOX,  l5HftECUft$I0N  PftOBS, 

iL_17HC0NDIT  TONAL  PROBS*__  15Xjr_  13HCUM  RECURSION*  12X,  *1.  -  CUHS, 

2  9X •  SMEAN  •TWEENS) 

00182  P  =  AlPHA*(Gb(l>  ♦  GB(1))  ♦  BETA  *  {  BG  (  1  >  ♦  BB(1)> 

00183 _ Q  «  P _ 

00184  ft  1  *  !♦  —  P 

00185  51  *  Q  *  FMUO  /  R1  - 

00186  PRINT  507*  P,  0*  Rl,  SI 

00187  507  FORMAT  (12X,$0*$,  8X,  E16.9,  1H*,  34X,  E18.12,  1H*,  9X, 

1  E 1 0  •  4  ,  S*S,  9X,  b 10  •  4  *  S*S) 

00188 _ P  BAR  =  0  - _ _ _ _ 

00189  PYAR  =  _ 

00V9G  DO  53  I-2,ft5 


00191  1S  _  «  I  -  1 

00192  P  =  ALPHAMGGU)  ♦  GB  (  I  )  )  ♦  8  E  T  A  *(BG(I>  ♦  BB  (  I  )  ) 

00193  Q  =  P  ♦  0 

00194 _ QH  =  1.  -  Q _ 

00195  R  =  P  /  R1 

00196  S  =  0  *  FMUO  /  QH 

00197  FIS  *  FLOATUS) 

00198  FIST*  FIS**2 

00199  P  BAR  =  PBAR  ♦  P  *  FIS 

00200  P VAR  =  P VAR  ♦  P  *  FIST 

00201  _  PRINT  501,  IS,  P,  Rf  6V  OH,  5 

66202  53  CONTINUE 

00203  501  FORHA?_t  8X*_fS,  \H*f  21  9x,  £15.8,  !«♦),  9X,  E18.12,  |H*,  9Xt 


1  E10.4,  1H*  ,  9X,  E10.4,  $*S> 

00204  606  CONTINUE 

00205  VAR  *  P  VAR  -  PBAR*PB AR 


00206 

66207 

00208 

00209 

00210 


06211 

00212 

00213 

00214 


SVAR 


SORT  <  VAR) 


CMEAN  =  PBAR  /  ftl 

CAU.  Tf*£ 


PRINT  305,  PBAR, 
305  FORMAT  (//  25X, 

1  25X 


TIMEX 


3000  CONTINUE _ 

if ( IK  .EQ. 
GO  TO  404 
ENO 


25x, 

25X, 


VAR,  SVAR,  CMEAN 
SPARTIAL  MEAN 
SAPPROXIMATE  VARIANCE 
$ APPROX l MATE  ST.  DEV. 


=  S, 


E 13 • 6/ 
El  3. 6/ 


1, 

SPARTIAL  CONDITIONAL  MEAN  *  S,  E13.6//1X,  A9//I 


0)  GO  TO  604 
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